Energetics of holes trapped in DNA 
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We report on our study of the electronic properties of guanine traps in the DNA surrounded by 
adenines. We have shown that for a typical range of DNA parameters, formation of the bound state 
of two holes at the same guanine trap is possible for the GGG and GGGG traps if the hole-hole 
interaction is weak, which can be achieved for the DNA in solutions. The origin of the two-hole 
bound state is the competition between the Coulomb repulsion and the phonon mediated attraction 
between the holes. For the hole-phonon coupling constant » 1 two holes will be at the same trap if 
the on-site hole-hole repulsion energy is < 0.9 eV. 

PACS numbers: 87.14.Gg,87.15.-v 



The motion of a hole on DNA is usually studied by 
detecting its trapping at a series of guanine sites incorpo- 
rated in the DNA Q. The hole can be injected into DNA 
through oxidation process. Oxidation occurs under UV 
irradiation and in the presence of certain oxidants. Since 
the guanine has the lowest oxidation potential 0, , the 
process takes place mainly at the guanine pair, which 
results in the formation of a guanine radical cation. A 
typical value for the oxidation potential difference Acm 
for the G-pair and the A-pair is about 0.1-0.5 eV 
which is also in the range of reported numerical estimates 
4]. Finally, after the oxidation process takes place we 
get a guanine radical cation, or a hole injected into DNA 
at the G-pair. The next step then is the migration of 
thepositive charge (hole) to the spots with lower energy 
[EE! which are the sequence of many guanines, e.g. GG 
or GGG. The importance of this migration is based on 
the fact that accumulation of holes at G-sequences makes 
these spots vulnerable to mutations. 

There are a few established pictures of hole migration 
through DNA until it reaches the spot with the lowest 
energy: (i) Direct hole tunneling between G and GG (or 
GGG) spots. This transport mechanism works only when 
the distance between the G and GG pairs is small (a few 
base pairs). The transport rate then depends exponen- 
tially on the transport distance 0,13- (ii) Hole migration 
through subsequent tunneling: This picture is applicable 
when the distance between the G and GG spots is large 
(about 10-15 base pairs) and when between the G spot, 
where the hole was originally injected, and the GG (or 
GGG) spots, there are a few G-pairs. In this case the 
hole migration is via subsequent tunneling between the 
G-pairs Q. Since this transport process is duffusive in 
nature, the transport rate has a power dependence on the 
distance between the G and GG spots, (iii) The polaronic 
hopping transport: This picture works when between the 
G and GG (GGG) spots there are only A-pairs. When 
the distance between G and GG is small the transport 
mechanism is again via tunneling (case (i)). However, 
for large distances, the hole at the G-spot tunnel into 
the nearest A base and then migrate via polaron hop- 
ping (diffusion) through the A. . .A sequence to the GG 



(or a GGG) spot pj. The transport rate has a power 
dependence on the distance between G and GG. 

Finally, the hole will be localized at the spots with the 
lowest energy, i.e. GG or GGG or GGGG (whichever 
is present in the DNA). A natural question that arises 
then is, if instead of a singe hole many holes are injected 
into the DNA, what will be the final configuration of 
the system? Which positions of the holes will have the 
lowest energy? Can the GG, GGG, or GGGG spots trap 
two holes simultaneously? In this case such spots will be 
highly reactive and vulnarable to mutation. Formation 
of the bound state of two holes trapped by G-sites is 
analogous to bipolaron formation in the homogeneous ID 
system 9J. In this paper, we report on the properties of 
the G-traps, surrounded by adenines. The main question 
we are addressing here is that under what condition the 
G-traps can accumulate multiple holes. 

The Hamiltonian of the hole system consists of three 
parts: (i) the tight-binding Hamiltonian which includes 
the hole hopping between the nearest base pairs and 
on-site energies of a hole, (ii) the hole-hole interaction 
Hamiltonian, and (iii) the Holstein's phonon Hamilto- 
nian with diagonal hole-phonon interaction |lfj| 

H = Ht + Hi + Hph, (1) 

where 

Ht = ^2 e i a \,a a i,° ~ 1 E [ a i> a i+i,a + h.c. , (2) 
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(3) 



where aj )(T is the annihilation operator of hole with spin 
a on site i, t\ is the on-site energy of hole, bi is the 
annihilation operator of phonon on site i, t is the hopping 
integral between the nearest base pairs (sites), lo is the 
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phonon frequency, \ is hole-phonon coupling constant, 
and n i a — aj CT a,- ia -. The Hamiltonian [Eq. l[TJl] without 
the T ^n/, was studied earlier for a homogeneous system in 
Ref. pj. 

In the tight-binding Hamiltonian [Eq. > we assume 
that the site i can be either a adenine or a guanine. We 
then take the energy of the hole on the A site as zero 
energy, i.e. ca = 0, and the energy of the hole on the G- 
site to be negative, eo — —Aga < 0. In the interaction 
Hamiltonian TLi, we take into account only the Hartree 
interaction between the holes. The first term in Eq. © 
describes the repulsion between the two holes with differ- 
ent spin. The holes can then occupy the same site. The 
second term in Eq. (JSJ) corresponds to two holes with the 
same spin. To get the basic idea about the typical range 
of interaction parameters when many holes are trapped 
by the G-sites, we introduce a single-parameter interac- 
tion potential of the form 



V, 



V(i-j) 2 + l' 



(5) 



where Vb is the on-site repulsion between the two holes. 
The origin of the form of Vij in Eq. is the finite 
spreading of the hole on-site states. This spreading is 
about the distance between the nearest base pairs. Al- 
though the actual dependence of interaction potential on 
the separation between the holes is more complicated 0] 
than Eq. (JjjJ , this difference is not important for our anal- 
ysis since only the on-site interaction plays the main role 
in formation of the bound state of two holes [13] • In the 
hole-phonon Hamiltonian Tt p h, we include only the opti- 
cal phonons [l4j with diagonal hole-phonon interaction, 
and do not take into account acoustic phonons which re- 
sults in non-diagonal hole-phonon interaction 0,0], i.e. 
modify the tunneling integral. In Eqs. ©-0 we also as- 
sumed that the hopping integral t, the phonon frequency 
uj, and the hole-phonon coupling constant x, do not de- 
pend on the type of the base pairs (A or G) . 

The form of the total Hamiltonian, Eqs. leads 
to four dimensionless parameters which characterize the 
system: The nonadiabaticity parameter |l7j 7 = fiio/t, 
the canonical hole-phonon coupling constant [l7| A = 
X 2 /(2fiw£), dimensionless hole-hole interaction strength 
Vo/t, and the dimensionless difference between on-site 
energies of G and A, Sqa — AgaA- 

We numerically determine the eigenfunctions and 
eigenvectors of the hole-phonon system by exactly di- 
agonalizing the Hamiltonian Eqs. Q-I@J for a finite size 
system consisting of six base pairs (sites). We also in- 
troduce limitations on the total number of phonons |l8j , 
^i n ph,i < N maxi where n p h^ is the number of phonons 
on site i. To compare the energy spectrum of the systems 
with different number of holes, we keep the maximum 
number of phonons per hole the same for all systems. For 
the two-hole system the maximum number of phonons is 
N m ax — 16 and for the one-hole system N max = 8. 

Our finite size system contains six sites which are orig- 
inally adenines. We then introduce the G-traps with 




FIG. 1: (a) The ground state energy of a single hole in a trap 
contaning Ng guanines for t — 0.2 eV and Aga = 0.3 eV 
as a function of the hole-phonon coupling constant, A. The 
numbers next to the lines are the number of guanines in the 
trap, (b) The energy Ei t i + £1,4 of two holes for Aga = 0.3 
eV as a function of A. The numbers next to the lines are the 
values of tunneling integral, t. 
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FIG. 2: Ground state energy of two holes in the trap contain- 
ing Ng = 3 guanines (a) and Ng = 4 guanines (b) as a func- 
tion of the inter- hole ineraction strength, Vb, for Aga = 0.1 
eV (dashed line) and Aga = 0.3 eV (solid line). 



different number Ng of guanines, G, GG, GGG, and 
GGGG, in the middle of the system. For example the sys- 
tem with two guanines is AAGGAA. For different traps 
we calculate the energy of the ground state of the sys- 
tems with one and two holes. Denoting the correspond- 
ing energies as Ei_m g (for the one-hole system with Nq 
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FIG. 3: The average number of holes (a) and the average 
number of phonons (b) for a two hole system in a GGGG 
trap are shown as a function of the base index. The tunneling 
integral is t — 0.3 eV and the hole-phonon coupling is A = 
1. Dots and triangles corresponds to inter-hole interaction 
strength Vo = 0.8 eV and 1.2 eV respectively. 



guanines) and l$2,Na (f° r the two- hole system with Nq 
guanines) we write the condition that the trap with Nq 
guanines will accomodate two holes as 



E2.N0 < Ei,n g + 



(6) 



The meaning of condition JHJl is as follows: If two holes 
are injected into the DNA then the lowest energy of 
the system corresponds to the case when the holes are 
trapped by the same trap with Nq guanines. The condi- 
tion |J?J| will determine the critical value of the hole-hole 
interaction strength, V cr . That means for Vq < V cr two 
holes will be trapped by the same trap with Nq guanines. 
For Vo > Vq t such a trapping is energetically unfavorable 
and two holes will be at different traps. 

As we shall see below the condition Eq. JSJ is not the 
condition for formation of bound state of two holes within 
the trap containing Nq guanines. The condition for for- 
mation of bound state of two holes should give different 
values for the critical hole- hole interaction strength, Vq. 
The values V cr and V b would coincide only for an in- 
finite homogeneous system. In our case the system is 
finite, which results in an inequality V cr < V b . The fi- 
nite system in our problem is actually the finite number 
of guanine sites in the traps. 

For our investigation of the system [Eqs. Ql-l@J], we 
consider the following typical DNA parameters: 0.1 eV 

< t < 0.3 eV, 0.1 eV < A GA < 0.5 0, EU, 0.05 eV 

< hw < 0.1 eV. For the dimensionless canonical hole- 
phonon coupling constant we have taken the value A = 
1. For this coupling constant, the size of the polaron 
is about 2-3 base pairs. Our calculations show that the 
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FIG. 4: Energies -Ei,i + £1,4 and -£2,4 of a two-hole system 
are shown as a function of hole-phonon coupling, A, by dashed 
and solid lines, respectively. Tunneling integral is t — 0.3 eV 
and Aga = 0.3 eV. 



critical value V cr is very small (V cr « 0.1 eV) when two 
holes have the same spin and they can not occupy the 
same site. This small value of V cr also illustrate the 
fact that the phonon mediated attraction between the 
holes are largest when the holes occupy the same site. 
Therefore, in what follows we shall consider only the case 
of two holes with opposite spin. 

In Fig. la the ground state energy of a single hole is 
plotted as a function of the hole-phonon coupling con- 
stant, A, for different types of traps. For Awl the dif- 
ference between the bound state of a hole in G and GG 
traps is about 0.03 eV, which is smaller than the value 
(0.05 eV) obtained in Ref. OS The size of the polaron in 
our calculations is 2-3 base pairs depending on the values 
of t and uj. As a next step we calculate the ground state 
energy of two holes, Ei^ G + £1,1, one of which is bound 
at a G trap and another one at a GGG or GGGG trap. 
As an example, in Fig. lb the ground state energy of two 
holes with the second hole bound by GGGG is plotted as 
a function of A for different values of the hopping integral, 
t. 

In Fig. 2 the ground state energy E2 : n g of two holes 
bound in a single trap is plotted for Nq — 3 (Fig. 2a) 
and for Nq — 4 (Fig. 2b) as a function of the hole-hole 
interaction strength for different values of A cm and t. 
Here we notice that at some critical value V b of the hole- 
hole interaction strength, there is a change of slope in the 
E2.n g (Vo) dependence. This critical value corresponds 
to the condition that the two holes are bound in the G- 
traps, forming a bipolaron. The illustration of this fact 
is given in Fig. 3. In Fig. 3a the average number of holes, 
(nh} — ( n i,a) + { n i.-a}, is shown as a function of the base 
pair index for a GGGG trap and two different values of 
the hole- hole interaction strength, V). It is clearly seen 
that for Vo — 0.8 eV < V b the two holes are almost at 
the same G sites, while at Vo = 1-2 eV > V b the holes are 
away from each other. The corresponding distribution of 
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TABLE I: Calculated values of Vq t for various values of the 
dimensionless DNA parameters 



t (eV) 


hw (eV) 


Aga (eV) 


V cr (eV) 


7 


$GA 


v cr /t 


0.1 


0.1 


0.1 


0.42 


1.00 


1.00 


4.2 


0.1 


0.1 


0.3 


0.45 


1.00 


3.00 


4.5 


0.1 


0.1 


0.5 


0.46 


1.00 


5.00 


4.6 


0.1 


0.05 


0.1 


0.36 


0.50 


1.00 


3.6 


0.2 


0.1 


0.1 


0.42 


0.50 


0.50 


3.25 


0.2 


0.1 


0.3 


0.59 


0.50 


1.50 


2.95 


0.2 


0.1 


0.5 


0.62 


0.50 


2.50 


3.1 


0.3 


0.1 


0.1 


0.82 


0.33 


0.33 


2.7 


0.3 


0.1 


0.3 


0.78 


0.33 


1.00 


2.6 


0.3 


0.1 


0.5 


0.82 


0.33 


1.67 


2.7 


0.3 


0.05 


0.3 


0.60 


0.17 


1.00 


1.97 


0.3 


0.05 


0.5 


0.62 


0.17 


1.67 


2.07 



the cLver£Lge number of phonons {ripfo }, is shown in Fig. 3b. 

Another critical value of Vq is introduced by equation 
Eq. (JSJ). The competition between E^i + £i,jv G and 
E2,n g is illustrated in Fig. 4. Comparing the energies 
Ex : Ng + ^1,1 and ^2,N G f° r A = 1 and different values 
of t, Aca, and u one can determine U cr . The result is 
summarized in Table 1 for the GGGG trap. The cor- 
responding results for the GGG trap is about 0.1 eV 
smaller values for Vq t . The dimensionless parameters, 
7, 8ga, and U cr /*, are also given in Table 1. From this 
data we can conclude that within the present range of 
parameters the dependence of Vq t on Aqa is weak, and 
Vq t ' jt depends mainly on 7. This dependence can be 
approximated by a linear function as 

U cr ps 3 7 i + 1.6* ps Zhw + 1.6*. (7) 

We see from these data that depending on the param- 



eters of DNA, the critical hole-hole interaction strength, 
Vy r can range from 0.4 eV to 0.8 eV. Numerical anal ysis 
of the electron correlations in different types of DNA [IJ 
shows that the hole-hole interaction strength is around 
0.9 eV for A-DNA and 1.5 eV for B-DNA. Additional sup- 
pression of the inter-hole interaction by a factor of ps 0.7 
|21j can occur for DNA in solution, when hole-hole inter- 
acton is screened by polar solvent molecules and mobile 
counterions. Under this condition trapping of two holes 
by GGG and GGGG traps would be possible. 

In conclusion, we have shown that for a typical range 
of DNA parameters, formation of bound state of two 
holes at the same guanine trap is possible for GGG and 
GGGG traps if the hole-hole interaction is sufficiently 
weak, which can be achieved for DNA in solutions. For- 
mation of the two- hole bound state results in double pos- 
itive charge of G-traps, which should modify their chem- 
ical properties and increase the reactivity of these traps. 
Experimental study of the two-hole states in G-traps 
should also give additional insight on the hole-phonon in- 
teraction strength. Formation of the bound state of two 
holes at the G-traps requires a strong hole-phonon in- 
teraction, which should overcome the hole-hole Coulomb 
repulsion. In our calculations the hole-phonon coupling 
constant was A = 1 which is larger than the experimen- 
tally reported A ps 0.2 in Ref. [22j. Hence, experimental 
observation of the two-hole bound state should give an 
additional estimate for the strength of hole-phonon in- 
teraction. 
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Research Chair Program and the Canadian Foundation 
for Innovation Grant. 
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